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Abstract 


In this paper we derive exact quantum Langevin equations for stochastic 
dynamics of large-scale inflation in de Sitter space. These quantum Langevin 
equations are the equivalent of the Wigner equation and are described by a 
system of stochastic differential equations. 

We present a formula for the calculation of the expectation value of a 
quantum operator whose Weyl symbol is a function of the large-scale infla¬ 
tion scalar field and its time derivative. The quantum expectation value is 
calculated as a mathematical expectation value over a stochastic process in 
an extended phase space, where the additional coordinate plays the role of a 
stochastic phase. 

The unique solution is obtained for the Cauchy problem for the Wigner 
equation for large-scale inflation. The stationary solution for the Wigner 
equation is found for an arbitrary potential. 

It is shown that the large-scale inflation scalar field in de Sitter space 
behaves as a quantum one-dimensional dissipative system, which supports the 
earlier results of Graziani and of Nakao, Nambu and Sasaki. But the analogy 
with a one-dimensional model of the quantum linearly damped anharmonic 
oscillator is not complete: the difference arises from the new time dependent 
commutation relation for the large-scale field and its time derivative. 

It is found that, for the large-scale inflation scalar field the large time 
asymptotics is equal to the ‘classical limit’. For the large time limit the quan¬ 
tum Langevin equations are just the classical stochastic Langevin equations 
(only the stationary state is defined by the quantum field theory). 
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I. INTRODUCTION 


While the quasi-classical picture of the inflationary universe scenario, based on a Fokker- 
Planck evolution equation for the probability distribution of the inflation held, is almost 
complete (see the basic papers H ). the essentially quantum-mechanical features of inha- 
tion are now a subject of investigation. 

We would like to note recent investigations into this problem carried out by Graziani 
0-0. Nakao, Nambu and Sasaki |^, and Hu, Paz and Zhang [|^. (The aim is not to 


present the full list of papers on this subject but to mention those which are closest to our 
considerations.) 

In the dynamics of an inhationary scalar held in a de Sitter background is investigated 


on the basis of the extended version of the stochastic approach proposed by Starobinsky [Q . 
In this approach, the scalar held operator is split into the long wavelength mode and the 
short wavelength mode. This split allows the reduction of the operator equation for the 
scalar held to Langevin equations of order ^/h. 

In the series of papers of Graziani, dealing with quantum probability distributions 

and the dynamics of the early universe, the approach is based on the Wigner function and 
its evolution equation, which is the Wigner equation. The author concentrats attention on 
the large-scale inhation ,where the spatial variable of the inhation scalar held is removed by 
averaging over a causal horizon volume. It is established that the quantum description of 
large-scale inhation in de Sitter space is equal to the quantum mechanics of a one-dimensional 
dissipative system. In it is shown how quantum Langevin equations can be derived to 
any order of ^/h when they correspond to the Wigner equation expanded in powers of h 
and truncated at some power of h. However the Wigner equation for large-scale inhation 
presented in |9HT^ is not accurate and we will return to this point in our paper. 

In paper dealing with quantum Brownian motion Hu, Paz and Zhang have derived an 
exact master equation for a quantum open system, which is an extension of earlier results 
obtained by Dekker [^], and Galdeira and Leggett |^. As was established by Vilenkin 


, Linde and Starobinsky [Q, the basic stochastic equation for the large-scale quantum 
evolution of inhation is similar to Brownian motion, that is the quantum Langevin equations 
of order ^/h. So the investigation of quantum Brownian motion is closely connected to further 
investigations of statistical and quantum ehects at the early universe. 

In this paper we restrict our consideration of the inhation scalar held to investigation 
of the dynamics of the large-scale or coarse-grained (> causal horizon) scalar held in a 
de Sitter background. Following in the main line Graziani [pHT^ we describe the evolution 
of the large-scale inhation by the Wigner equation. Then the aim of this work is to derive 
exact quantum Langevin equations (to all orders of y/h), which describe the stochastic 
dynamics of the large-scale inhation. 

Thus, our approach is opposite in some sense to the approach of Hu, Paz and Zhang 
These authors move toward the possible observation of macroscopic ehects from the 
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search for an adequate description for statistical and quantum ehects, while our way is to 
start from the macro-level to obtain an equivalent stochastic description. 
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II. BASIC FORMULATIONS 


The Lagrangian density of the inflation scalar held ^{x,t) in a de Sitter backgronnd is 


C 






where the backgronnd metric is assnmed to have the form 

= —dt^ + o?{t)dx^, 


( 1 ) 

( 2 ) 


g denotes a metric, t is a time, a; is a three-dimensional spatial coordinate, V(<h) is a potential, 
and a{t) is a scale factor (some positive fnnction). 

Then the eqnation for the classical inhation scalar held ^{x,t) is 


d(t) d 
.dt"^ a{t) dt 



A, 


^x,t)+vmx,t)] 


0 , 


(3) 


where A^, is the spatial Laplace operator , the dot over a symbol means a time derivative 
y{t) = dy/dt, and V'[<h] = dV[<h]/d<h. 

The expansion of the nniverse is assnmed. 


d{t) 

a{t) 


H{t), 


(4) 


where H{t) is some non-negative integrable function of time. In other words, the scale factor 

is 

a(t) = a{0) exp{J iL(r)dr}, (5) 

where a(0) is the value of the scale factor at t = 0, the beginning of inhation. In this model 
the beginning of inhation is taken to coincide with the origin of the universe. 

Strictly speaking, in our consideration we deal with an expanding Friedmann-Robertson- 
Walker (FRW) universe, ‘approximatly’ close to a de Sitter universe in the sense: H{t) ~ 
Const. 

From (1) one obtains the canonical momentum conjugate to the held ^{x,t) '■ 


n(a;, t) = 


and the Hamiltonian density : 


dC 


d{d^/dt) 


= a(t) 


dt 


( 9 $ 


= la{t)-^U^ + ia(t)[V,<I>]2 + a(t)3V($). 


( 6 ) 


(7) 


When interested in large-scale (> causal horizon) physics, a coarse-graining procedure is 
utilised and it leads to a coarse-grained or averaged scalar held ^x(t) 
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( 8 ) 


V jQx 

where the index X is a label referring to the center of a region over which $(a;, f) is 
averaged, and V is its volume. The volume of spatial averaging is taken not smaller than a 
causal horizon volume: 


V" > (9) 

with the causal horizon (or the ‘coordinate horizon’ in terms of Ellis 0) given by 

i{t) = [ a(r)“^(ir. 

Jo 

For de Sitter space 


a(t) = exp{iff}, H = Const 


( 10 ) 


then i{t) = and the volume of averaging is chosen to be E = {4:/3)xH~^ . 

For ‘approximatly’ de Sitter space we will think about a volume of averaging as constant 
for all time. This always can be assumed if i{t) < Const. The techniques we use allow 
consider time dependent volume V{t) along the same lines but we shall not consider this 
case for the sake of brevity. (See remarks in Section EH) 

Each large-scale region (labeled by X) can be considered as a separate quantum- 
mechanical system because each lies outside of its neighbours’ light cone: there is no exchange 
of information between large-scale regions. The proht of the coarse-graining procedure is 
that it reduces the quantum held problem to a quantum mechanical problem. At the same 
time there are still some peculiarities, following from the held theory, which do not make 
the analogy with quantum mechanics complete. We will point out these peculiarities in this 
section. 

After performing the coarse-graining (^ the spatial varying term in the Euler- Lagrange 
eqation (^) can be neglected due to a smaller factor a~‘^{t)V~'^^^ (see [^) and the equation 
for the large-scale inhation held ‘hx(t) becomes 


^x{t) + 3H{t)^x{t) + U'[<!>x{t)] = 0, (11) 

where U{^) is the coarse-grained potential. 

The averaging of the momentum (j^) gives 

Ux{t) = a{tf^x{t). (12) 

The Lagrangian (density) for the coarse-grained held <hx(^) is 

£|4v]=a(«)='lii.5..-(/(4x)]. (13) 


This gives the Euler-Lagrange equation ( 0 ) and keeps the averaged momentum (|T^ canon¬ 
ically conjugate to <hx(t): 


Hxit) 


dC[^x] 

d^x 


(14) 
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Now the Hamiltonian (density) for the coarse-grained field ^x{t) is 


H(tx,Tlx; t) = + a(t)3t7(<*x), 


( 15 ) 


it can be considered as the classical Hamiltonian for the large-scale inflation held, in the 
sense that the equations of motion produced by this Hamiltonian 


d^x _ dH dllx _ dU 

~dr ~ ~ 


( 16 ) 


are equivalent to the held equation ([Tl|). 

To make the step from the classical equation (O) to the quantum equation it is necessary 
to quantize ( pd^ taking into account and ( 0 ) . This leads at the end to the quantum 
mechanics of a one-dimensional dissipative system in the description on <f)x, variables. 

One can apply the canonical procedure, which is based on quantal noise operators and 
conserves the fundamental commutator for canonical position and momentum operators in 
the course of time, or the inhuence-functional method of Feynman and Vernon |^. Both 
ways give the master equation for the ‘reduced’ density operator p (see |]l5| , pT| for techniques), 
describing the time-evolution for the large-scale inhation. 


dp 

m 


^[H,p|-.^o(«)®C(()|4x,|4>.v.p]]. 


(17) 


Here 7i is an operator form for the classical Hamiltonian ( 0 ). 1.1 stands for a commutator, 
and the dihusion coefficient D{t) is, in general, some non-negative function. In particular, the 
dihusion coefficient can be assumed to be D{t) = 3H{t)- Const, where Const is determined 
by physical parameters of the system at equilibrium. 

Some one who would follow, formally, the quantisation procedure in quantum mechanics 
would notice immediately two points by which the master eqation (O) differs from the 
formally obtained one: hrstly, the new constant h/V instead of Planck’s constant h, and 
secondly, the factor a{t)^ in the diffusion term, which reflects the explicit time dependence 
of the Hamiltonian (|T^. It is time to discuss how these peculiarities come from the quantum 
held theory. 

The origin of the scaling 


h h/V (18) 

in equation (p!7[), where V is a volume of averaging (^, is that the fundamental commutator 
for canonical position and momentum operators in the relativistic quantum held theory 

[$(a;,f),fl(?/,t)] = ihS{x-y) 

is transformed by coarse-graining procedure (P) to 

[^x{t),flx{t)]=th/V, (19) 

[«l>x(t),ny(t)] =0, V. 
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An other way to derive the master equation (pAj) would be to use the following momentum 
and Hamiltonian: 


Unit) = VUxit), 


( 20 ) 


Wn = VH(<i>x,nx-,t). 


For momentum Hn the fundamental commutator is 

[^x{t),Un{t)]=th. ( 21 ) 

The description in terms of 0-0 is natural for the large-scale region Returning back 
to the spatially dependent held the Hamiltonian for H is found to be 


Hn= [ Uix,t)dx, 
Jn 

or, in view of the Hamiltonian (density) 


nn = vni^x,Ux;t) 


( 22 ) 


(23) 


(in (^31) the spatially varying term of the original held [Vx*h]^ was neglected ). 

Now let us discuss the dihusion term in the master equation (|^). It should be mentioned 
that the general form of the dihusion term in the master equation is 

where C{t) is some time-dependent coefficient. This coefficient is dehned with respect to 
a suitably chosen vacuum state in the held theory. For de Sitter space it is the so-called 
Bunch- Davies vacuum [Q- In terms of ([T7|) this means that it is determined by the 


physical parameters of the system at equilibrium. 

If in accordance with the held theory (of inhation) we assume that the energy density of 
the equilibrium state is an invariant in de Sitter space then 


{■Hn(t))/aiffV = {H{i>x,tlx;t))/a(tf = Const, 


(24) 


where, to obtain the energy density, the Hamiltonian Tin is divided by the proper volume 
ait)^V. The expectation values, ( ) , of the operators are taken on the stationary solution of 
the master equation. For ‘approximatly’ de Sitter space. Hit) is neglected so the condition 
(1^) can be assumed. 

To fulhl (P^ we have obtained the factor a(f)® in the dihusion term of (P!7|). Taking into 
account that the constant in (^) is proportional to ih/V), in (^) we can show explicitly 
the dependence 

C(t) = ^a(t)‘D(,t), 

which is valid for large-scale inhation. We will use master equation (0 as a starting point 
for investigation of the stochastic dynamics of the large-scale inhation. 

From this point onwards in the paper the large-scale scalar held is refered to as *h(t), 
omitting the index X. 
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III. THE WIGNER EQUATION FOR LARGE-SCALE INFLATION 


The Wigner function W(q, p; t) p^ , |23[| is a function on the classical phase space and de¬ 
scribes the distribution of position and momentum. The Wigner function is not a probability 
distribution since it can assume negative values; W{q,p;t) is a real function. 

As a density matrix the Wigner function contains all of the information corresponding 

to the quantum state. The expectation value for an arbitrary operator y4(<h,<h) can be 
calculated by the formula 


(i) = / dqdpA(q,p)W{q,p;t), 


(25) 


where A{q,p) is a Weyl symbol for the operator 24(<f), d>). 

Let us derive the time-evolution equation for the Wigner function for large-scale inflation, 


that is equivalent to the master equation (pT]) . 

For the Wigner operator W the equation has the same form as for the density operator 
P (O) (because coefficients a{t) and D{t) depend only on t) 


dW U - - U « 

~m = ■ Ti “(*> !«• 

where Weyl symbols for W and Ti are 

W ^ W(q,n;t), 

n ^ n{q, n) = + a{tfv{q). 


(26) 


(27) 


Here we assume the following correspondence between phase space variables (g, H) or {q,p) 
and variables of the large-scale inflation scalar held (<h, <F): 

<h = g, = p, a{tf^ = n. (28) 

To obtain an equation for the Wigner function W{q, H; t) from its operator form (|2B|) we 
will use the formula for the composition of operators in the Weyl formalism. Due to this 
formula and the commutator relation ( [TP|) a Weyl symbol A{q, H) for the composition 
of two operators A = BC is dehned via the correspondence 




^ d II 

,n 


^,11 ih ^ d II d , 

= 


(29) 


where B{q, H) and C{q, H) are Weyl symbols for the operators 13 and C, and the numbers 
I,II over operators show the order in which the operators act. 

Using (|27|) ,(|2^ ) one can obtain the following correspondence between operators and their 
Weyl symbols: 


^ ^ ? /? f) 

[n,W] ^ -ya{t)-^U—W{q,U-t) 


+a{ty 


U q + 


ih d 


-U q 


ih d \ 


W{q,U;t), 


(30) 
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[■j.B.w]] 


(31) 


Combining relations (PUD-(^) in accordance with equation (P^ one has the time evolu¬ 
tion equation for the Wigner function W{q, If; t) 


dW(q,Il;t) .dW h , ,d'^W 


V / 33 


U[q + 


ih d \ 

wan) 


Uiq 


ih d \ 

wm) 


W(q,-a-,t) 


(32) 


Equation (^) is called the Wigner equation. (Note that the commonly used way to derive 
the Wigner equation is to apply the coordinate representation for the Wigner function 


1 ft ft 

W{g,n-,t) = -I e‘”-(q-—x\f,\q + —x)dx 


(33) 


to the master equation (p^Tl). Actually it is a longer way because it deals with multiple 
integrations and integrations by parts while the Weyl formalism gives an answer in a plain 
algebraic way.) 

Because we are interested in the distribution of <h and <h let us make the transition from 
W{q,Il;t) to W{q,p;t) by relations (1^) and equality 

(34) 


W(q,p-,t) = a(t)^W(q,Il-,t) 


which follows from (|25|) . For the Wigner function W{q,p;t) one has the time evolution 
equation 




^ 2 33 

where H(t) is defined by ( 
The Wigner equation 


dt 




-rr, / 3-3 9 


W{q,p;t), (35) 


can be rewritten in the following equivalent form: 




+ Tpa{t)^ duW{q,p - u;t) ■I{q,u;t), 

zft J—00 


where 


1 /■°° 

= — dye 

Ztt J—00 


—luy 


2V 


+ NT?®W 7rr7«(^) y 


2V 


(36) 


(37) 


Equation (P5|), or equation (|56D, is the complete Wigner equation (to all orders of h) for 
the large-scale inflation in de Sitter space. It describes the time evolution for the distribution 
of <h and <f) in a sense of equality (p^). The expansion of the Universe, which is described 





















by the scale factor a{t) and equality (H), gives the origin for a dissipation term, with the 
disspation coefficient 3H{t), in the Wigner equation. In this sense large-scale inflation can 
be considered as a quantum one-dimensional dissipative system, which supports Graziani’s 
statement |Q. At the same time the expansion of the Universe influences the potential terms 
of equation (^) (or (|3^)), which was missed in 

In [|I^|l4| the Wigner equation for a quantum disspative system was derived for a har¬ 
monic oscillator’s potential U{q) = a;^g^/2. For an arbitrary potential U{q) a truncated 
Wigner equation, or a Kramers-Moyal equation, is commonly used in the literature instead 
of the complete Wigner equation (see for example 


TB[]). In our case (|55|) it would be 


dW{q,p;t) 

dt 


dW , ,,,, h 

-P^r~ + ^H{t) — {pW) + 


dq 
dU{q) 


dW 


dq dp 


dp 

+ 0{h^ 


V 


dp'^ 


(38) 


where 0 {h^) is a value of order h^. 

The potential term in equation (|35|) can be expanded in powers of h by Taylor’s series, 
which leads to equation 


E 

A;=l,3,5,.. 


k\ \2V^ ) dq^ dp^ 


(39) 


In the series of papers [PHT^ the expansion in powers of h of the Wigner equation was 
used to show how to improve the accuracy of (P5|). 

The potential term in Wigner equation (|35|) (or in the equivalent representations (^)- 
(P7|) and (PP)) on a macro-level contains information about quantum noise on the micro-level, 
when 


d^U{q) 

dq^ 


^ 0 for fc > 3, 


(40) 


This quantum noise is non-Gaussian noise. Gondition (|40|) is always fulhlled when a potential 
deviates from a harmonic oscillator’s potential. The diffusion term, containing d‘^W/dp‘^, 
in the Wigner equation also represents quantum noise (Gaussian noise) on a macro-level. 
Some authors consider a possibility that the Wigner equation may contain second order 
derivatives d'^W/dq^, d^W/dqdp, which are also of quantum origin (These terms 

would correspond to additional terms [If, [If, p]], [*h, [If, p]], [If, ["h, p]] in the master equation 
(17).) These terms have no physical sense because they lead to wrong Langevin equations: 
dqt/dt 7^ pt or d^/dt ^ 4). 

In this paper we will work with the complete Wigner equation (|36|) -(|37D for large-scale 
inflation to derive equations for the corresponding stochastic dynamics. Let us supply the 
evolution equation (|36D with the initial condition 


lU(g,p; 0) = IFo(g,p), 


(41) 
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where the initial Wigner function is chosen to satisfy the following properties: 

a) hho(<?)P) is a real function, 

b) JWo{q,p)dqdp=l, 

c) J[Wo{q,p)]^dqdp < ^a{0f, 

which follows from general properties for the Wigner function (see |p3|| ). 


(42) 

(43) 

(44) 


IV. QUANTUM LANGEVIN EQUATIONS 

In this section we deduce quantum Langevin equations which are equal to the complete 
Wigner equation but we start with some assumption about the potential U{q). 


UJ 


U{q) = + 

where p{p) is a bounded measure such that 


exp{-iqp} ij(dp), 


(45) 


J p{dp) < Const, 
p^ti{dp) < Const 


(46) 


The hrst term in the right hand side of equation (^5]) is just the harmonic oscillator’s 
potential while the second term can be considered as a deviation from it. The potential is 
assumed to be real. 

To include more model potentials used in the theory of inflation o;^ is allowed to be either 
positive or negative, or to be zero. A parabolic potential connected with ‘chaotic inflation’ 
is included in (^) with p{dp) = 0. A double-well potential, representable in the form (^Sf ) 
is, for example 

U{q) = uj^q‘^/2 + K cos{kq) ■ I[-{ 37 r/ 2 k),( 3 n/ 2 k)] (q), (47) 

where uj,K,k are some parameters (real values) and /a( o') is the identihcator of a set A. 


Uiq) = 


1, liqe A 
0, if g ^ A. 


Form (^) rules out model potentials which are polynomials higher than second order in 
q. Remember here that model potentials of polynomial type have appeared in the theory of 
inflation through expansion of C(q) in powers of q near q = 0 (see, for example, |^). Thus, 
almost all physical potentials for the inflation held <h(t) can be represented as 
For potential (^) the Wigner equation becomes 


dW(q,p;t) 
dt 


= -P- 


dW 

dq 


+ UJ q- 


dW 


d 


h 


dp +mt)-{pW) + -D{t) 


d^W 


iV P ^3 
a(t)3 


V dp^ 

n 


dp{p') exp{-iqp'} w(^q,p+ —a{t) V; 


h 


-W{q,p- ^a{t) ‘^p'-t 


(48) 
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QiPiP' ^ ^ ^ [0)^]) supplied by the initial condition 

W{q,p;0) = Woiq,p), 


( 49 ) 


such that properties (|4^)-(^^ are satished. 

In order to solve this equation and to deduce stochastic equations, which describe a path 
in a phase space {q,p) = ($, $), we should reduce it to the form of a forward Kolmogorov 
equation 0( p.l02). Such an approach was hrst proposed by Chebotarov and Maslov |2^ 
and developed much further by Combe, Guerra et ah 0- 

To reduce equation (^8]) to a forward Kolmogorov equation we need to transform the 


last term in the right hand side of it to a standard form. It can be done by introduction of 
a new function 


r 21/ 1 

f{q-,P-,0-,t) = W{q,p-,t)expl —^ ^ R^/mod27r, 

where 6* is a new variable, and function A{t) will be dehned later. 

Inserting (|^) into (^8]) one has 


(50) 


df{q,p,0]t) df 2 9f , ^ 


dt 


= -Pvp + w q-?^ + 3H{t) — {pf) + TrDit)- 

oq op op V Op^ 


I + ^p'u,e + ^u- qp;t^ 


21 / 

X 


■ [5(m + 1) + S{u - 1)] 

2 VdA{t) 


-fiq^p^^A) 

r 

-a{tff{q,p,GA) ■ J pW) - 


h dt 


S(a,p, 0\t), 


( 51 ) 


where 5{u) is a Dirac delta-function. 
Let us put 


A{t) = J a{T)^dT ■ J p{dp). (52) 

Note that for scale factor (|I^) 

A{t) = J p{dp) ■ (e^-^* - 1)/3R. (53) 

Let mn{dpdu]t) be the following measure: 

V n 

mn{dpdu;t) = —a{ty^p{dp)[6{u -l- 1) -f- S{u — l)]du. (54) 

The measure mri{dpdu;t) is a bounded measure on (D (—1) ® (-1-1) multiplied by time- 
dependent function a{t)^. For an arbitrary function 'ijj{p,u) an integral over this measure 
is 

I i,(p,u)mn{dpdu-,t) = ya(tfl[i,{p,+l)+i,(p,-l)]p(dp). ( 55 ) 
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Now, with (^) and (Q), equation (^Tj) becomes 


df{q,p,0;t) 

dt 




dp"^ 


-fiq,p,0;t) 


mn{dp'du] t), 


( 56 ) 


which is the forward Kolmogorov equation. 

If equation is formally supplied by the initial condition 

f{q,P,0;t) = 5{q - qo) s{p - Po) 6{9 - Oo), (57) 

then, according to the theory of stochastic differential equations (SDE), there exists a three- 
dimensional stochastic process 

6 = (<!>*, <i>i,0i) (58) 

for which the function f{q,p,0-,t) is a probability distribution. This means that for an 
arbitrary function h{q,p, 6), which is continuous and periodic with the period 2n on variable 
6 and for hxed 9 belongs to the class C^’^(R^ 0 R^), 


J Kq^P,d)f{q,p,9]t)dqdpd9 = Eh(^t), 
where the symbol E denotes mathematical expectation and 

it=Q = (q'o,Po, 9o) 


(59) 


(60) 


(see equality (^)). 

Using the generalized Ito formula for stochastic differentials (see [^,p.270, formula (13)) 
one deduces an SDE for the stochastic process @) 


d^t = ^tdt, 

= — [3H (t)$t -l- -|- 

-^a(t)"^y puun{dpdu-,dt), 




1/2 


dwf 


dQt = 


71 


u + <htp 


Vfi^dpdu] dt). 


(61) 


(62) 

(63) 


Here Wt is a one-dimensional Wiener process (or Brownian motion) and un{dpdu]dt) is a 
Poisson measure on R' (g) (—1) ® (+1) ® [Oj T] non-homogeneous with respect to translation 
on [0, T] such that 


Eil^nidp du; dt)] = mn{dpdu;t)dt. (64) 

Qt is an additional stochastic variable, which can be interpreted as a stochastic phase (see 
equation (|50D ). It appears only due to deviation of the potential (]45|) from the harmonic 
oscillator’s potential. 
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In order for the stochastic process (|^ - (|63D to have a unique solution on a time interval 
[0, T], right continuous with probability 1, it is enough to have conditions (|4^) , and to assume 
that functions a{t), H{t), D{t) and their hrst derivatives are continuous functions. 

Let us denote the stochastic process (|^-(|6^) with initial condition (|57|) by 

^t{qo,Po,do) = {^t{qo,Po,do),^t{qo,Po,do),Qt{qo,Po,do)). (65) 

If the initial condition for equation (^61) 

f{q,p,0;0) = fo{q,p,e) (66) 

belongs to the class of generalized functions then instead of ( |5D|) one has 

J h{q,p,e)f{q,p,e;t)dqdpde = 

= J 'Eh{^t{q,P,d),^t{q,P,d),Qt{q,P,d))fo{q,p,9)dqdpde. (67) 

Now one can readily hnd the correspondence between an integral over the Wigner func¬ 
tion W{q,p-,t), governed by equation (|48|) , and the stochastic process (|6l|)- (|6^ . Let us 
assume that 


fo{q,P,d) = Wo{q,p) 

and let the function h{q,p, 6) in ( |57|) be of the form 

h{q,P,9) = h{q,p)exp{-ie}. 
Then from (^^ and (|50|) one hnds 


h {q,p)W{q,p;t)dqdp = 

2V 


= exp 


{^7l(f)} JE{h{<l>t{q,p, 0), ^t{q,P, 0)) exp[-iQt{q,P, 0)]}Wo{q,p)dqdp, (68) 


where J d9 = 27i was used and the function A{t) is dehned by (^^. 

Formula (|6^ gives one the expectation value of a quantum operator A(<F, <h) with its 
Weyl symbol h{q,p) (compare ( |68[) with (p5|)). 

Equations (|^ ) - (|6^) are exact quantum Langevin equations for the large-scale inflation 
in de Sitter space associated with the master equation (P!7|). 

In 113 an attempt was made to derive the quantum Langevin equations for potentials 
with a polynomial growth higher than second order in q. As was mentioned, such potentials 
are excluded in our consideration, which is restricted by (^). In the deduction of the 


quantum Langevin equations is based on a general expansion in powers of h of the Wigner 
equation, as in our representation (|3^. From []^ it follows that for polynomial potentials, 
of order higher than second, the quantum Langevin equations can be derived exactly only 
to order \Jh/V (corresponding to the truncated Wigner equation of order fi/V). Already 
the hrst correction of order h/V to the quantum Langevin equations cannot be calculated 
precisely (there is no explicit representation for the noise terms in the quantum Langevin 
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equations). The origin of this problem is related to the term SP’W/dp^ in the truncated 
Wigner equation of order {Ti/VY. 

Let us rewrite the quantum Langevin equations (|^)-(|63|) in an integral form to show 
explicitly dependence on initial data (g,p, 0): 


^t{q,P,0) = q+ ^r{q,P,0)dT, 


^t(g,P,0)=p- / [3H{T)^riq,P,0) + u^^riq,P,0)]dT + 
Jo 

h 


2 hlV 2 rt 


V 


D(T)^^^dvJr 


fl 

~wrP / a{T)~^p'uun{dp'du-,dT), 

2V Jo 

0 f(?,p,O)^|^ J ^u + ^r{q,p,0)p' //fi(rfp'dt£;dT)|/mod27r. 


(69) 


(70) 

(71) 


Equations ®-(0) describe a stochastic path in phase space (<h, <h) starting at 

a point {q,p) when t = 0. The stochastic phase 0* plays role in the hnal formula (^) for 
the expectation value of a quantum operator and can be interpreted as a contribution of the 
stochastic path <hi(g,p),$(g,p). 

Equation 


together with equations (|^) -(|6^ can be used for numerical simulations 

to calculate the expectation value for quantum operator A(<|), <h), which is much simpler 
than solving the Wigner equation because SDE (|^ - (|63D are of hrst order. 


V. SOLUTION OF THE WIGNER EQUATION 


The aim of this section is to show how the Wigner function itself can be expressed by the 
expectation value with respect to a stochastic process in the extended phase space {q,p, 0). 
In such an expression is found for the case 


dW 1 .. .. 


(72) 


while in the case under consideration (^),(]^) the diffusion term and the expansion of the 
Universe play an essential role. 

To hnd a solution for the Wigner equation (^8|) one should transform it to a backward 
Kolmogorov equation |^( p.300). This can be done by introduction of the following func¬ 
tion: 

f{q,P,d]t) = W{q,p;t)expl^-3j^ iL(r)dr - (73) 

where 9 G R^/mod27r, and A(t) is the same function (^) as in Section B 
The time evolution of / is governed now by the equation 


df{q,P,9;t) 

dt 


d‘^f 

= -P^ + + 77^ Witt 


dq 


'5 


h 

V' 


dpi 


+ J f(q,P + ^a{'t) yu,e -^u- qp']t 


-f{q,p,d;t) 


mn{dp'du] t), 


(74) 
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which would be a backward Kolmogorov equation for the backward time t', by setting 
t = T — f, where 0 < t < T. The measure mf,{dpdu; t) is dehned by (|5^) . 

If equation ( 0 ) is supplied by an initial condition 

fiq,p,9;0) = fo{q,p,0), (75) 

where /o is a continuous function, then the solution for O can be represented as 

f{q,P,9-,t) = Efo{QT{q,p,0-,T -t),PT{q,p,e-,T -t),QT{q,P,d]T -t)), 0<t<T, (76) 


where the stochastic process can be found by applying the generalized Ito formula for stochas¬ 
tic differentials (p.270), 


Qs{q, p,0;T -t) = q- j Priq, p,0;T - t)dT, 

JT-t 

Ps{q,p,9]T -t) =p+ f [3i7(r)P^(g, p,9;T -t)+ uj^Qr{q,P,9;T - t)]dT 

JT-t 


(77) 


+ 


2hih2 PS 


V 

h 


iT-t 


D^rY^'^dWr 


+7717 / f ^p'uun{dp'du;dT), (78) 

zv JT-tJ 

Qs{q,p,9]T -t) = {9 - ^u + Qr{q,p,9;T -t)Y Uh{dp du] dr)^/mod2n, (79) 


where 0<T — t<s<T and unidp'du] dr) is a Poisson measure on ® (—1) ® (-1-1) ® [0, T] 
with the intensity dehned by (|64|). 

Returning to the Wigner function one has 

W{q,p-,t) = exp^^S J i7(r)dr-7 ^A(f)| ■ 

■E{lKo(QT(g,P, 0; T - f), PT(g,P, 0; T - t)) exp[i0T(g,p, 0; T - t)]}, (80) 

where the stochastic process is dehned by equations (|77D-(|79D with 9 = 0. 

Conditions for existence and uniqueness of the solution for SDE ([77|) - (|79D are the same 
as for SDE (|69D-([7TD. 

Let us assume additionally that 


Jp^p{dp) < Const, /c = 3,4, 


(81) 


and that the initial Wigner function WQ{q,p) is twice continuously diherentiable in p and 
once in x, and that its hrst and second order partials are bounded. Then the Wigner 
function (^) is twice continuously diherentiable in p and once in x, diherentiable in t and 
is the unique solution for the Cauchy problem ([7^)-(|75[). 

Let us consider how properties (^ ) -(|4^ ) for the initial Wigner function are preserved in 
the course of time. 

a) The reality of the Wigner function is preserved. To prove this proposition it is nec¬ 
essary to use the decomposition of the Wigner function into a diherence of two positive 
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functions (see formula (4.4)). The group of transformation (|80|) preserves this decom¬ 
position. 

b) The normalization for the Wigner function 


J W{q,p-,t)dqdp = 1, t > 0, 


is fulhlled due to the equality 


d 

dt 


W{q,p-, t)dqdp = 0 


which follows directly from the Wigner equation. 

c) The presence of the dissipation and diffusion terms in the Wigner equation destroys 
the inequality (H )- If now becomes 


W‘^{q,p;t)dqdp < J WQ{q,p)dqdp ■ J if(r)(ir| < ^^a(t)^. (82) 


To derive (|82|) it is necessary to take a time derivative of the expression on the left hand 
side of (|8^) and to use the Wigner equation. After this one has 


W^{q,p; t)dqdp = 


a{t) 


a(0)J 

2h A 

'V Jo 


J Wo{q,p)dqdp 

a(0)i3 


a(r) 


D{t) 


•r dW{q,p-,T] 
V dp 


dqdp 


dr}. 


(83) 


and (|8^) follows from 

If there is no diffusion, D(t) =0, for such a system, starting at t = 0 from a pure state 


[[Wo{q,p)]'^dqdp = 7r^a(0)^ (< 7w^a(0)^ for a mixed state), 
J zirh znh 

it is possible to follow the pure state in the course of time due to equality 


[ W‘^{q,p-,t)dq dp = (< 7w^a(t)^ for a mixed state). 

J zirri zTih 


(84) 


However, diffusion {D{t) ^ 0) smears the picture and one cannot disinguish pure and mixed 
states by inequality (|8^) . 

For the Wigner function W{q,Il;t) (see (P^)) relations (P3|), (|5^) become 


W^{q,U;t)dqdU= W^{q,lV)dq dU 


V 


W^{q, U;t)dqdU < 


dW{q,U;T) 

an 


dq dU 


dr, 


2Tih 


(85) 

( 86 ) 
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VI. LARGE TIME ASYMPTOTICS 


If H{t) ^ Const the scale factor a{t) increases exponentially with time and for large time 
the Wigner equation becomes 




(87) 


to order a{t) what follows from 

The expectation value at large time T for an operator A with Weyl symbol h{q,p) is 


{A)t = J h{q,p)W{q,p-,T)dqdp 

= jEh{^Tiq,P, 0), <i>T(g,P, 0)) Wo{q,p)dqdp, 


where the stochastic process (<l)t, 4)^) is governed now by 


d^t 

d^t 


^tdt, 


-[3H{t)^t + U{^A]dt + 


-2h 

.V 


1/2 


B(t) 


dwt 


( 88 ) 


(89) 


with the initial condition (<l)o, <ho) = 

Note that the stochastic differentials ([89| ) are equivalent to the quantum Langevin equa¬ 
tions (0-® for the scalar held <h(t) in the large time limit (not for the beginning of the 
inhation if the potential U{q) satisfhes the unequality (^0|) . 


VII. STATIONARY STATES FOR LARGE-SCALE INFLATION 


Stationary or equilibrium states for large-scale inhation are described by the stationary 
Wigner function W{q,p) = limt^oo governed by the Wigner equation 


dW dUiq) dW h 

~ ^ dq ' Up ^ V°^*'^lpF 


(90) 


for an arbitrary potential U{q). This equation has time dependent coefficients coming from 
the time dependence of the scale factor. 

For the stationary Wigner function W{q,p) one can easily hnd the following expectation 
values for the operators: 

^ {i) = 0, 

{(/(4)) = o, 

From the last equality one can see that the dihusion coefficient D{t) should be 

D{t) = 3H{t)a (94) 


(91) 

(92) 

(93) 
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(at least for large time), where the constant a is defined by the choice of stationary state: 


( 6 ) 



The stationary Wigner function is found to be 


W(<i,p) = Afexp|-(y + (7 ((/))v'/(R(t)|, 


(95) 


(96) 


where N is the normalization constant. 

For a particular case, with the scale factor given by (pTl) and potential U{q) = uj‘^q‘^/2, 
the Bunch-Davies vacuum is given by 


3H^n UH 
^ ^ ~ 8 ^ ~ 

which yields a = H/2 'k and D = 3H‘^/2n. 

For this case the stationary solution is 


= ^exp{-(p2+o;V)^}, V 


dvr 


(97) 


(98) 


where one must assume that 


H >uj/tt 


for the inequality (^df ) to be satisfied by the Wigner function (p8|). 


(99) 


VIII. CONCLUSIONS 


1. The appearance of the dissipation term 3H{t){d/dp){pW{q,p;t)) in the Wigner equa¬ 
tion, after transition from phase space ($, fl) to (<F, <h), supports the earlier result of Graziani 
^ and Nakao, Nambu and Sasaki |^, that the large-scale inflation scalar held behaves as 
a quantum one-dimensional dissipative system. Nevertheless, this analogy is not complete: 
it is destroyed by a new commutation relation 

[^,h] = ya{t)-\ ( 100 ) 

where a{t) is the scale factor in de Sitter metric and it rehects the expansion of the Universe 


Comparing the Wigner equation (|48|) for the large-scale inhation (or more general case 
(^ ) - (l^) ) with the Wigner equation for a quantum one-dimensional linearly damped un¬ 
harmonic oscillator (see EHig). one can see that the expansion of the Universe amplihes 
the role of the potential term, which is a deviation from the harmonic oscillator potential. 
Explicitly this amplihcation appears in formula (|68D for the expectation value of a quantum 
operator and in formula (|80|) for the Wigner function, where function A{t) (|5^) has the 


factor J^a^T^dr (without expasion this factor would be 1). At the same time , as t —>■ oo, 
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the commutation relation (|101|) leads to degeneration of the jump process, in the quantum 
Langevin equations, into a continuous process (see Section P. 

2. As a consequence of our investigation we have the following statement: for the large- 
scale inflation scalar held the asymptotic t > cxo is equal to the classical limit. 

In the limit h ^ 0 the Wigner equation (or (|i8|)) turns out to be 


at oq op oq op 


( 101 ) 


and the corresponding Langevin equations are 


d^t = ^tdt, 

d^t = -[SH{t)^t + U'{^t)]dt. (102) 

Equations ( |102|) are just the classical deterministic equations of motion, equivalent to the 
held equation (pT]). 

However, ( |101| ) is not the ‘classical limit’ for the Wigner equation (^)- (^) (only the 
truncated Wigner equation of order ho. In the classical limit the dihusion term in the Wigner 
equation (H)-® is not proportional to h/V because, instead of (P3|) for a stationary state, 
one has 




(103) 


The ‘classical limit’ of the Wigner equation for the large-scale inhation scalar held is 


Ot Oq Op Oq Op Op^ 

with Dciit) = 3H(t)aci. The classical stochastic Langevin equations are 


(104) 


d^t = ^tdt, 

d^t = + U\^t)]dt + [20,1(1)]^/“^ dwf (105) 


One can see that that in the ‘classical limit’ the Wigner equation and Langevin equations 
(( |104j) and ( |105|) ) are the same as those of the large time asymptotic ( equations (|87|) 
and (^)). For complete coincidence the coefficient of the second derivative in the Wigner 
equations (|57| ) and ( |1U4|) should be presented in the form 3H{t){^‘^)st, where the expectation 
value is taken on the corresponding stationary state. 

3. In this paper each large-scale region (for the coarse-graining procedure) is considered 
as an independent quantum mechanical system. If, nevertheless, it is necessary to take into 
account interaction with the environment, then for linear interaction the master equation 
for the ‘reduced’ density operator p is 


dp 

m 


+^(l<i'./5ni - |n,/5 4>l) - ja{t)0D{t)l^, P.p]], 


(106) 
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where X{t) is the dissipation coefficient (originated by interaction with the environment), 'H 
is an operator form for the Hamiltonian ([I5| ) of the system without dissipation, and {, } 
stands for an anti-commutator. 

In phase space ($, <h) the Wigner equation corresponding to (p.U6|) is 


dW{q,p-,t) 

dt 


dW 



d‘^W 



+ |3/f(«) + Mt)]-^{pW) + XD(t) X 

1/ ^OO 

—a(f)^ / duW{q,p — u;t)X{q,u;t), 
iTl J — OO 


(107) 


where X(g, M;f) is dehned by (pTf) . 

Equation (|107| ) is an equation of the same type as the Wigner equation (pbl). Thus, (|107|) 
can be treated as the Wigner equation considered in this paper. 

4. Since in our consideration on phase space (<h, 4)) the scaling ha{t)~^/V depends already 
on t, it is not a problem to take into account a time dependence of the coarse-graining volume 
V{t). The scaling becomes /V(t), where H(0) = Const ^ 0 is assumed. 

Therefor the result can be extended on an expanding FRW space-time. The principal 
point in this case is the new time dependence of the canonical momentum conjugate to 


nn(«) = V(t)a(t)^ix(t) 


(108) 


and of the Hamiltonian: 

nn = IV(t)-^a(t)-^n^^ + V(t)a(t)^[/(<^x) (109) 

(compare with (pOD and (|^) ). Also the condition (|^) , taken for ‘approximatly’ de Sitter 
space, should be replaced by a condition corresponding to each concrete FRW space-time. 
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